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Abstract. Observations show that about the 20% of the Universe is composed by 
invisible (dark) matter (DM), for which many candidates have been proposed. 
In particular, the anomalous behavior of rotational curves of galaxies (i.e. the 
flattening at large distance instead of the Keplerian fall) requires that this mat- 
ter is distributed in an extended halo around the galaxy. In order to reproduce 
this matter density profiles in Newtonian gravity and in cold dark matter (CDM) 
paradigm (in which the DM particles are collisionless) , many ad- hoc approxi- 
mations are required. The flattening of rotational curves can be explained by a 
suitable modification of gravitational force in bigravity theories, together with 
mirror matter model that predicts the existence of a dark sector in which DM 
has the same physical properties of visible matter. As an additional result, the 
Newton constant is different at distances much less and much greater than 20 kpc. 
[To be published in Proc. Int. Conf. "Atomic Clocks and Fundamental Constants" 
ACFC 2007 - Bad Honnef, Cermany, June 2007] 



1 Introduction 

Cosmological observations show the Universe to be nearly flat, i.e. the energy density is very 
close to the critical one: f2tot — 1- We can separate the different contributions present in the 
total density. Mainly, we can write Htot = + ^A, in which we distinguish the contributions 
due to matter {{2m — 0.24±0.02) and to the dark energy (cosmological term, = 0.76±0.02). 
In particular f2M = ^b + ^d, where Qb = 0.042±0.005 and Qn = 0.20±0.02 are, respectively, 
the baryonic and dark matter components 1 . 

Let us focus on the dark matter problem. CMD model is in agreement with experimental 
results on Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) that confirm 
the presence of the dark matter component. Many DM candidates have been proposed in the 
literature as axions (m ~ 10~^ eV), neutralinos (m ^ TeV), wimpzillas (m ^ 10^^ GeV) and so 
on. However, the question why the fractions of the baryonic and dark components are so close, 
flo/^B ~ 5, remains unresolved. 

A most convincing proof for DM comes from rotational curves of galaxies and cluster dy- 
namics. In the Newtonian picture, at large distance from the center of a galaxy one expects 
that the velocity along a circular orbits behaves as v{r) cx l/^/r (also known as Keplerian fall). 
This is due to the fact that the gravitational potential of a galaxy outside the inner core (bulge) 
is ~ GM/r, where G is the Newton constant. Instead, one observes that in these regions v 
is approximately constant. In order to explain this anomalous behavior, without modifying the 
Newtonian paradigm, one has to suppose the existence of dark matter distributed in a extended 
spherically symmetric halo around the galaxy, according to the "isothermal" mass distribution 
profile p{r) oc (1 + (r/a)^)^^ , where a is a scale radius. Since gravity is universal between visible 
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and dark matter, a point-like source of both types of matter generates a potential 

^(r) ^ -{Ml + M2) (1) 
r 

where Mi and M2 are, respectively, the visible and dark components. CDM model assumes 
that the DM particles are coUisionless. Due to this property, N-body numerical simulations tend 
actually to predict a different behavior: the DM distribution has a cusp profile p{r) cx 
with a = 1 1.5 }2|3| . For many galaxies these cusp profiles do not reproduce the observed 
rotational curves as well as the isothermal profile. 

An alternative possible candidate is the Mirror matter [4]. According to this model our 
Universe is made of two similar gauge sectors. In other words, in parallel to our sector of the 
ordinary particles (0-) and interactions described by the Standard Model, there exists a hidden 
sector (M-) that is an exact duplicate of ordinary sector in which particles and interactions have 
exactly the same characteristics, and the two sectors are connected by the common gravity 
(see for review [51). Therefore, if the mirror sector exists, then the Universe should contain 
along with the ordinary particles (electrons, nucleons, photons, etc.) also their mirror partners 
with exactly the same masses (mirror electrons, mirror nucleons, mirror photons, etc.). Mirror 
matter, invisible in terms of ordinary photons, can naturally constitute dark matter. One should 
stress that the fact that O- and M-sectors have the same microphysics, does not imply that 
their cosmological evolutions should be the same too. In fact, if mirror particles had the same 
temperature in the early universe as ordinary ones, this would be in immediate conflict with 
Big Bang Nucleosynthesis (BBN). The BBN limit on the effective number of extra neutrinos 
implies that the temperature of the mirror sector T' must be at least about 2 times smaller 
than the temperature T of the ordinary sector, which makes mirror baryons viable candidate for 
dark matter. In particular, the mirror dark matter scenario would give the same pattern of LSS 
and CMB as the standatd CDM if T'/T < 0.2 or so [B]. In addition, the baryon asymmetry 
of the Universe can be generated via out-of-equilibrium B — L and CP violating processes 
between ordinary and mirror particles [7] whose mechanism could explain the intriguing puzzle 
of the correspondence between the visible and dark matter fractions in the Universe, naturally 
predicting the ratio Hd/^b ^ 1 10 [8J. 

However, in contrast to the coUisionless CDM, mirror baryons obviously constitute coUisional 
and dissipative dark matter. Therefore, one should expect that mirror matter undergoes a 
dissipative collapse and thus in the galaxies it is distributed in a similar manner as the visible 
matter instead of producing extended quasi-spherical CDM halos. Indeed even the hidden sector 
undergoes a dissipative collapse as the visible sector that follows an exponential profile p(r) oc 
g-i/fo^ In this way the distribution of the dark matter is more compact in the center of the 
galaxy and is not extended as the CDM halo. 

Since gravity is universal between the two sectors, this mirror dark matter hypothesis gets 
into difficulties to explain the flat rotational curves of galaxies. However we can suppose that 
each sector has its own gravity and that mixing term produces a suitable modification of gravity 
at large distanceQ. In particular we show that the interaction term allows us to obtain a massive 
graviton and leads to a modified potential. A test mass of type 1 at distance r from the origin in 
which there is a sources of both types of matter {Mi and A^2), instead of ([T]), feels a potential 

0(r) = ^(A^i+A^2) + ^V^(Xi-X2), (2) 
Ir Zr 

where G is the Newton constant and is the range of the massive graviton. Notice that at 
small distance r the test mass interacts only with Mi through the ordinary Newton 

potential, whereas, at large distance r ^ the test mass interact with the sum of the two 
kinds of matter {Mi + M-i)!"!. This result, together with the mirror matter hypothesis, enables 
us to reproduce the observed rotational curves of galaxies. 

^ The particle mixing phenomena between ordinary and mirror sectors were discussed in the literature 
for photons [9], neutrinos [10], neutrons [11], etc., as well as possible common gauge interactions between 
two sectors [12) . The mixing between the ordinary and mirror gravitons was first discussed in our recent 
papers [19I2U| . 
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2 The Model 



Let us consider a theory with two dynamical metrics (71,2 /ii/, each of them interacting with its 
own matter. The total action contains two Hilbert-Einstein actions and a mixed term V: 



,Ml ^ , ^,M? 



VFii^Ri + >Ci) + V9~i{^R2 + C2) - m'(5i.92)^/'V(5i, 52) 



2 v.^^x 2 



(3) 



where Afi_2 are the Planck masses (in general different) and £i_2 are the corresponding matter 
Lagrangians (respectively, ordinary matter or type 1 and dark matter or type 2) . The action ([2|) 
describes a more generic bigravity theory, that is the simplest case of the multigravity theory 
which considers N metrics interacting each other through a mixing term (see for review [13j). 
The interaction term breaks down the invariance under the diffeomorphism group Di (8) D2 
to a diagonal diffeomorphism -Di+2- The two metrics 171 and 172, in the flat (Minkowski) space 
approximation, can be written as 51,2/ii/ — 'ni,2p,u + ^i,2/ii//Afi,2- The mixed term in ([3]) can 
induce the non diagonal rank 1 mass matrix between two gravitons hi and /12 which has one 
massless and one massive eigenstates: 



hf^u = cos-dhifj^iy + sini?/i2^i/ 



(4) 



where 1? is a mixing angle: tan?? = M2/M1. In the case Mi = Af2, i-e. 'd — 7r/4, the rotation 
(|4]) reduces to even and odd combinations of hi_2fj.u- The massless state ft,^^ is the ordinary 
graviton that exhibits a Newtonian potential ^ l/r universally coupled with both matters. 
The massive state h^i, in turn can have a Lorentz breaking (LB) mass pattern [14] 

-Cmass = ("^0^00 + '^'n^l^Oi ~ ^l^ij + ^l^i " ^mlJioohi^ , (5) 

(0 and i = 1,2,3 are the time and space indices, respectively) that can induce, for a suitable 
combination of the masses m^'s {14, a Yukawa term in the potential ^ (l/r)e~'"/''"* in weak 
field limit approximation. 

From these solutions, in general it follows that a test particle of ordinary matter (type 1) 
feels a static potential induced by a point-like source containing the mass fractions Mi and 
M.2 of the two types of matter, as 

0W = :r (-^1 + -^2) + (tan^ §Mi - M2) , (6) 

Zr Zr 

where G ^ l/(87r(Aff + M|)) and ^ and Vm are parameters that depend on the pattern of the 
masses in the ([S]). The first term is mediated by massless gravity and the second term by the 
massive one. The symmetric case, i.e. Mi = M2, allow us to simplify the potential ^ in the 
following form 

GMi fl + ^e-'-^'-^-\ , GM2 fl-^e-''^'- , 

that shows directly the modification with respect to the Newtonian law ^ . Let us distinguish 
two cases: the mass term ([5]) is Lorentz invariant or Lorentz breaking. In the first case the only 
consistent mass term (without ghosts) is Pauli-Fierz type [15] (i.e. mo = and mi, 2, 3, 4 — rn). 
In this case in ^ one has 1^ 4/3 and therefore a deviation from the Genral Relativity pre- 
diction for the light bending in gravitational field 0, also noted as van Dan-Veltmann-Zakharov 
discontinuity (vDVZ). 



^ In this case one can define the discontinuity parameter: 5 = 1 + (sin^ i^)/3. Experimental limits on 
the post-Newtonian gravity [17] requires 5 = 1.0000 ± 0.0001, therefore -ff ~ 10"'^. As a consequence 
the Lorentz invariant case Mi = M2 is excluded. 
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In the second case, in general the masses in ([5]) are different. The discontinuity is absent if, 
in the mass term, toq = 0, mi ^ and/or m2 ^ m^. However, the bigravity ^ produces a 
mass term with mi — and therefore does not give a Yukawa potential [18] . The addition of 
a third auxiliary metric gs, which works as a bridge between gi and g2 and decouples with a 
Planck mass M3 3> Mi^2, allows us to have an effective bigravity theory with mi ^ [TOj. In 
this case the massive gravity introduces a Yukawa potential (with range l/{\^m4)) without 
vDVZ discontinuity (i.e. ^ = 1) as in ([2]). This potential, in the mirror matter model, can 
explain the rotational curves of galaxies with the condition A42 — WMi and r„i ~ 20 kpc (i.e. 

Ho/nM ^ 10). 

3 Rotational Curves of the Galaxies 

The galaxy rotational curves describe the velocity of stars and interstellar gases as a function of 
the distance r from the center. For the sake of simplicity we apply our model to disk galaxies, 
where most of the matter (about 2/3) is concentrated in the inner region, called bulge. Indeed, 
one can suppose that the matter density along the profile follows approximately the luminosity, 
which decreases exponentially with r moving out from the center. In a spherically symmetric 
approximation the visible matter has an exponential distribution depending only on r, i.e. 
p(r) = poe^^^^°, where po is the density in the central region and ro is the size of the bulge. 
The Newtonian theory, that takes into account only visible baryonic mass, does not explain 
the flattening of the rotational curves at large distance and in general needs to introduce an 
extended halo composed by dark matter with different density profile with respect to the visible 
one. 

In bigravity theory, in which we consider the mirror matter as dark matter candidate, a 
different explanation emerges ([20]). Let us assume the following hypotheses 

— The ordinary and dark matters interact only via gravity, modified according to the bigravity 



— The two types of matter have similar density profiles in the galaxy, i.e. exponential along 
the disk, p{r) oc exp(— r/rg). 

The velocity of an object at distance r is determined by equating the centrifugal acceleration 
with the radial component of gravitational acceleration a(r) derived from the potential ([2|). For 
instance, from the gravitational field of point-like source of both types sitting in the center, we 



In order to obtain the total force on a star moving approximately along a circular orbit around 
the galaxy center we have to integrate ([8|) on the matter spatial distribution. Notice that for 
r << r^, the influence of matter of type 2 on a particle 1 is negligible and the behavior is 
Newtonian. The behavior in the opposite limit r » Vm is also essentially Newtonian, though 
the test particle feels the presence of the total mass Mi + M2 and the effective Newton constant 
is G/2. In the region r Vm there is a significant deviation from the Newtonian theory due to 
the presence of matter of type 2, resulting in a enhancement of v. This result avoids the cusp 
problem which is present in the context of the CDM paradigm and reproduces a isothermal-like 
shape {Md/-Mb ~ 5 and pbm{i") oc (l-|-(r/ro)^)~^). For instance, let us consider agalaxy with 
Ml = IO^^Mq and the bulge size tq ~ 3 kpc (e.g. the Milky Way). In Fig. [1] we compare the 
rotational curves fitted with the isothermal DM halo (a ~ 8 kpc) in standard gravity and those 
fitted with the exponential DM profile in the bigravity theory (ro = 5.4 kpc for the invisible 
distribution). Notice that both curves have approximately the same behavior. In addition, we 
show the visible matter contribution to the velocity reproducing the Keplerian fall proportional 
to 

The flat rotational curves can be reproduced varying the parameter Vm and the mass ratio 
M2/M1. Fig. [2]shows the rotational curves for different M2/M1, with rm = 20 kpc. Notice, 
that for Ml = M2 we obtain the Keplerian fall as we expect from the Newtonian potential ([2]). 
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Fig. 1. Rotational curves fitted by the isothermal DM halos in standard gravity (dashed) and by the 
exponential DM profile in the bigravity theory (solid) with M2I M\ ~ 10 and r™, = 20. Contribution 
of the visible matter to the velocity given by the dotted curve. 
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Fig. 2. The rotational curves for different M2/M1, with r™, — 20 kpc. Notice, that for A4i — M2 we 
obtain the Keplerian fall. 



Incidentally, if both the sectors have the same content of matter the potential ([2]) is indistin- 
guishable from the standard Newton potential </) = G{Mi/r) generated by the visible matter. 



4 Conclusions 

If the Universe is made of two separate sectors, one visible and one hidden, each of them can 
have its own gravity and the two metrics can interact as in a bigravity pattern. This lead to 
a large distance modification of the gravitational force, though it remains Newtonian at small 
distances where each type of matter feels only itself, in agreement with the precision test of 
General Relativity in the Solar System. 

The mirror matter as dark matter, together with Lorentz breaking bigravity can explain the 
flattening of rotational curves of galaxy at large distance. This model supposes that both dark 
matter and modification of gravity are present. The main advantage of this theory is that, due 
to the same property of both types of matter, the dark matter can have a mass distribution in 
galaxy similar to the visible one. This avoids the need to invoke the presence of extended halo 
distributions that are in conflict with numerical simulations in the CMD paradigm. 

As an additional result, it interesting to note that the Newton constants for type 1 - type 1 
attraction differs by a factor 2 for large and small distances: G]y{r <C r^) — G and GAr(r ^ 
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Tm) = G/2. In the general potential ([6]) the Newton constant can be measured as GN{r ^ 
rm) — G[(l + ^tan^i!?)/2] between type 1 - type 1 matter at small distance, while at large 
distance we have GN{r ^ r„i) — G. Moreover, it is interesting to note that, if ^ > 1, at small 
distance r ^ one has antigravity between type 1 and type 2 objects. 
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